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We calculate the geometric phase of a spin-1/2 particle coupled to an external environment com- 
prising N spin-1/2 particle in the framework of open quantum systems. We analyze the decoherence 
factor and the deviation of the geometric phase under a nonunitary evolution from the one gained 
under an unitary one. We show the dependence upon the system's and bath's parameter and ana- 
lyze the range of validity of the perturbative approximation. Finally, we discuss the implications of 
our results. 
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I. INTRODUCTION 

A quantum system undergoing a cyclic evolution behaves quite differently from its classical analogue: 
the state vector retains the history of its evolution in the form of a geometric quantum phase. Berry [l[ 
demonstrated that closed quantum systems could acquire phases that are geometric in nature. He showed 
that, besides the usual dynamical phase, an additional phase related to the geometry of the space state is 
generated during an adiabatic evolution. A phase only dependent upon the area covered by the motion of 
the system. 

A renewed interest in geometric phases (GP) has recently blossomed by the proposal of using GP for 
quantum computation. The use of fault tolerant quantum gates in quantum information requires the study 
of the geometric phase in open quantum systems 0, H, 0, M, H 0, Si- This has been motivated by the fact 
that all realistic quantum systems are coupled to their surrounding environments. No matter how weak the 
coupling that prevents the system from being isolated, the evolution of an open quantum system is plagued by 
nonunitary features like decoherence and dissipation. Decoherence, in particular, is a quantum effect whereby 
the system looses its ability to exhibit coherent behaviour and appears as soon as the two interfering partial 
waves shift the environment into states orthogonal to each other Q. Nowadays, decoherence stands as a 
serious obstacle in quantum information processing. 

Starting with the seminal paper of Zurek [lo| , several authors analyzed the decoherence due to a collection 
of independent spins. In this paper we discuss a two-level system coupled to a one-dimensional array of spin- 
1/2 particles. This model evolves reigned by a Hamiltonian which encompasses Ising and Universality in one 
simple approach. In some trivial cases, we can solve the problem exactly. Most of the works done so far are 
based on the so called "central spin model" , where the two level system is coupled isotropically to all spins 
of the bath, which implies a enormously simplification of the derivation. To ensure decoherence, the central 
spin (the system) is coupled to the transverse component of the bath. 

Working towards having a realistic description of geometric phases, we have introduced field decoherence 
by coupling the central two level system to an external environment. In this paper we investigate the 
geometric phase of a spin-1/2 particle interacting with a spin environment. The model is the most simple 
case of an Ising chain and a further consideration to the Zurek model, since we are considering the self 
energy of the environment spin-1/2 particles. We calculate and analyze the decoherence factor induced 
by the coupling to an external bath. Within the general framework of GPs in open quantum systems, we 
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estimate the corresponding GP for the spin- 1/2 particle in the weak coupling regime. We further compare 
this result to the general GP obtained numerically and analyze its dependence upon the bath and system 
parameters. 

This paper is organized as follows. In Sec. II, we present the model and analyze the decoherence factor 
due to the interaction with the environment. In Sec. III. wc estimate the geometric phase for the central 
spin and thoroughly study its behaviour by performing numerical and analytical calculations. Finally, we 
conclude our results in Sec. IV. 



II. MODEL 



In this section we examine the decoherence induced by disordered interacting spin baths at finite tem- 
perature. Our choice of the bath is the most simple case for an Ising chain so as to facilitate an analytical 
study of the decoherence and the geometric phase for the model. A further study of the geometric phase for 
an Ising chain with random spin-spin interactions will be studied elsewhere. Our system is comprised by a 
central spin CTc coupled to a bath of N Ising spins through the interaction Hamiltonian 



N 

J7i = a,^®^A.af. (1) 



The coupling between the central spin and the bath is characterized by the coupling constant A. The presence 
of a non interacting spin chain will be modeled by a bath Hamiltonian such as 

N 

HB = J2oj,al (2) 

1=1 

where erf denotes the Pauli operators of the bath spins. The self Hamiltonian of the central spin is the 
common one Hg = h/2Qa^. This model results a further approach to the model studied in (lo| since we are 
considering not only the interaction Hamiltonian but also the sclf-Hamiltonians of the N spins of the bath. 
In addition, we must note that [Hb, Hi] ^ 0. 

Wc start with a total separable state of the form such as 

|M'(0)) = (a|0)+/3| l))x Ixo), (3) 

where | 0) and | 1) are the state eigenstates in the ct^ basis and xo is the initial state of the bath. At a later 
time, the total state of the system is 

|vl/(t))=e-'^*|*(0)), (4) 

where H is the evolution total Hamiltonian. The main object in the study of decoherence is the reduced 
density matrix which is obtained after tracing out the environment's degrees of freedom. In our case, the 
reduced density matrix may be written as 

where F(t) = TrB[e-^(^«+^)VB(0)e*(^«"^)*]. We have named V = X^^i ^i^t after the interaction Hamil- 
tonian has operated on the central spin. 

We can choose a pure state for the initial state of the bath as |xo) = Otiil'^ilOi) + Allj))- Since 
Pb(0) = |xo)(xo|, the decoherence factor F{t) yields 



^(*) = n (l - sin\^oj^ + X^t)) . (6) 
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We can see that this decoherence factor does not depend on the external temperature. It is also possible to 
check that F{0) = 1. As it was done in llOl - we can evaluate the mean value of the decoherence factor. This 
expression is similar to that obtained in but for a different initial state of the environment. As in that case, 
{F{t))T^oo 0, when T ^ oo. If we estimate the average dispersion as = {F{t)'^) — {F{t))^ = J^iLiPi' 
with Pi = (1 — ^^^a^^a^ )■ Under the assumption that all pi are approximately equal, the average fluctuations 
from zero are 




Therefore, large environments can effectively induced decoherence on the central spin system. 



III. GEOMETRIC PHASE ENVIRONMENTALLY CORRECTED 



In this section, we shall compute the geometric phase for the central spin and analyze its deviation from 
the unitary geometric phase for a spin-1/2 particle. 

A proper generalization of the geometric phase for unitary evolution to that for non unitary evolution 
is crucial for practical implementations of geometric quantum computation. In [T^ . a quantum kinematic 
approach was proposed and the geometric phase (GP) for a mixed state under nonunitary evolution has been 
defined by Tong et.al. as 

$ = arg{^ v/£fe(0)£fc(r)(ffc(0)|vl',.(r)) x /o^ (8) 

k 

where ek{t) are the eigenvalues and l^ffc) the eigenstates of the reduced density matrix p,- (obtained after 
tracing over the reservoir degrees of freedom) . In the last definition, r denotes a time after the total system 
completes a cyclic evolution when it is isolated from the environment. Taking the effect of the environment 
into account, the system no longer undergoes a cyclic evolution. However, we will consider a quasi cyclic 
path V : t e [0,t] with r = 27r/i7 (O the system's frequency). When the system is open, the original GP, 
i.e. the one that would have been obtained if the system had been closed is modified. This means, in a 
general case, the phase is $ = + S^, where 5$ depends on the kind of environment coupled to the main 
system 0, [ij, [l^. For a spin-1/2 particle in SU{2), the unitary GP is known to be <I>^ = 7r(l + cos(6'o)). 
It is worth noticing that the proposed GP is gauge invariant and leads to the well known results when the 
evolution is unitary. 

To this end, we shall start by choosing a pure state of the form: 

|*(0)) =cos(0o/2)|O)+sin(0o/2)|l). 

This is the same to have assumed a ~ cos(6'o/2) and (3 = sin(6'o/2) in our above deduction of the decoherence 
factor. Then, for times t > 0, the state of the system is 



|vl/+(i)) =e-f"cos(0+(t))|O)+sin(0+(i))|l), (9) 

where we have defined tan(6'+(t)) = ta.n{9o/2)F^^{t). It is easy to check that, when F{t) = 1, we re-obtain 
the results for the unitary case. From Eq.Q, one can easily write the reduced density matrix Pr{t) and 
estimate its eigenvalues and eigenvectors, as was done in Ref.0]. In order to estimate the geometric phase, 
we only need the eigenvector |^+(t)) since £-(0) = 0, and, hence, the only contribution to the phase comes 
from that eigenvector and its corresponding eigenvalue (see Eq.®). In the case we are considering here, and 
after performing the time derivatives of Eq. we can note that the eigenvalues and eigenvectors arc real 
functions. Consequently, the geometric phase is 



$ = O /" cos^{0+{t)). 
Jo 



(10) 
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As we can deduce from the relation between 9+{t) and the decoherence factor, Eq. fTU]) is not easily 
computed. In Fig.l, we show the behaviour of the GP as a function of the initial position on the Bloch 
Sphere (^o) smd the value of the coupling constant (A) for an external environment comprising 10 spin-1/2 
particles. 




FIG. 1: Exact GP as a function of the coupling constant and the initial angle in the Bloch sphere. For simplicity, we 
have assumed that all the bath spins have the same frequency and coupling constant and A'^ = 10. 

In order to achieve an analytical result, we shall assume all the bath spins to have the same coupling 
constant and frequency. Then, the decoherence factor becomes F{t) = YiiLi fi{t) — f{t)^ a-nd we can forget 
about the product function. Besides, we shall perform a perturbative expansion in powers of the coupling 
constant A. In such a case, 

$ = 7r(l + cos(6lo)) + ivf-) sin^{eo)(TT - -^smiATTLu/n) 



u) J V 4a; 



(11) 



Notice that the lowest correction in \/lo of Eg. pTj) is quadratic, which means that in case of low de- 
coherence we recover the unitary GP up to the first order in This reflects the resilient of the phase 
against the environment. If the decoherence rate is sufficiently small, the probability amplitude of staying in 
the same eigcnstate is almost stationary. This amounts to effectively project the state back to the original 
eigenstate, whenever it tends to be driven away by the environment. Therefore, the state trajectory on the 
projective Hilbert space tends to be unaffected by the decoherence (up to first order), thereby leaving the 
area enclosed by the path, and hence the geometric phase, unchanged. Similarly, this robustness against 
decoherence can be observed in different two- level systems 0, [^, Q . 

In Figs m and [31 we show the comparison between the exact calculation fEq. pU)) ) of the GP and the 
perturbative approximation fEg. pT]) ) for different values of the coupling constant. In Figl^l we plot the 
nonunitary GP for a case of a small environment consisting of = 10 spins. The lined curves, whether solid 
or dashed, correspond to exact calculations of the GP, for different initial position on the Bloch sphere. The 
perturbative estimations are plotted with lines and points in each case. Therein, it is possible to see that for 
a small initial angle, i.e. small value of 6q, there is almost no significant difference between both approaches. 
However, for an initial state of bigger value of the environmentally-induced correction to the unitary GP 
becomes important which is shown in Figl3l for an environment of = 100 spins. This difference among 
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FIG. 2: Comparison between the exact GP and the perturbative approximation for different initial positions on the 
Bloch sphere and A*' = 10 spin environments. The lines (whether solid or dashed) represent the exact GP while the 
lines with different shaped points are the corresponding perturbative perturbative calculations. 




FIG. 3: Comparison between the exact GP and the perturbative approximation for different initial positions on the 
Bloch sphere and N = 100 spin environments. The lines (whether solid or dashed) represent the exact GP while the 
lines with different shaped points are the corresponding perturbative perturbative calculations. 



the approaches becomes more relevant as the number of spin environments increases and can be noted by 
comparing the scales of the y-axis in each figure. This can easily be imderstood by looking at Eg. pTj) . 
Therein, one can identify a new "effective" coupling constant X^g = iV(A/w)^, which means that one must 
consider the product of both quantities and that the perturbative approximation will hold for a compromise 
between them. That is to say that for larger environments (bigger values of A^), the perturbative estimation 
will remain valid for smaller values of the coupling constant. However, this estimation remains valid for a 
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bigger range of the coupling constant if we reduce the number of spins comprising the environment. 

Another feature that can be studied is the deviation from the unitary value of GP for both different 
approaches: exact and pcrturbativc, as shown in FigUl Therein, the GP vs the initial position on the Bloch 
Sphere (^o) is plotted for different environment sizes and coupling constants. The lines (dotted or solid) 
correspond to the exact calculation of the GP while the lines with different shaped dots are the corresponding 
perturbativc calculations for each case. As can be expected, when the coupling between the system and the 
environment increases, the perturbative calculation of the GP differs from the exact one. This difference is 
more relevant for larger environments, as can be observed in FigUthe comparison between an environment 
comprising N = 10 spin-1/2 particles and one with N = 100 (and the same coupling constant A = 0.05). 
Not only, is it easy to note that this deviation is bigger for larger environments but also, it can be observed 
that it is more relevant for smaller angles as the number of environment spins increases. 
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FIG. 4: Deviation of the exact GP and the perturbative approximation from the unitary value of the GP for different 
initial positions on the Bloch sphere and for different number of spin environments (TV — 10 and A'^ = 100) and 
different values of A. 



Finally, we can analyze if the phase retains its geometric nature. Even though we know that the evolution 
of the system is not suppused to be cyclic due to the presence of the environment, under certain external 
conditions it may result a valid assumption. To that end, we plot the exact GP of Eq.(l5]) as a function of 
the number of cycles the system evolves normalized by the exact GP of only one cycle. Therein, we can 
note a nice topological feature of the phase: the phase is dependent on the number of times the path is 
traversed, i.e. the winding number. In the end, the GP adquired by the system will be the same. However, 
this does not mean that dephasing will not affect the meassurement of the GP. The visibility fringes of any 
meassurement made on the system will be reduced due to the presence of the environment [3, Q . 



IV. CONCLUSIONS 



So far we have studied the corrections to the geometric phase under a non unitary evolution induced by 
the presence of decohcrence. The model is a further consideration in the spin-spin dccohcrcncc model since 
we are also considering the self energy of the spin-1/2 particles comprising the environment. In addition, it 
is the simplest approach to study the nonunitary geometric phase in an Ising chain. 

As expected the geometric phase is affected by decoherence. This correction is quadratic in the coupling 
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FIG. 5: Exact geometric phase of m cycles normalized to the exact geometric phase of one cycle as a function of the 
number of cycles performed by the system. The plot results independent of Acff . 



constant, so this affection is negligible when we work in the weak coupling constant X/ui ^ 1. However, the 
correction is linear in the environment's size A'^ as well, which might have a considerable contribution for 
larger environments. 

We have also shown an analytical expression for the environmentally induced correction to the GP, based 
in a series expansion in powers of the system-environment strength. What is worth noticing is that due to 
the dependence with the initial angle ^Oi for small angles, there is not a significant correction induced by the 
environment, almost all the geometric character of the phase is in the unitary term. On the contrary, for a 
bigger value of the corrections to the unitary term become important. This conclusion is stressed as long 
as the number of spins in the environment increases. On the other hand, we have analyzed the validity of the 
perturbative calculation of the induced correction to the unitary term, as it is shown in Fig[4l Therein, we 
have shown that the larger environment, the shorter the values of A for which the perturbative expansions 
hold. 

Finally, we have studied the geometric nature of the total GP (the one that includes the environmental 
correction) changing the single cyclic evolution ruled by r = 27r/0 to m times r, with m an integer number 
(the winding number). We have shown (see FiglS]) that the total GP can be seen as m times the GP for 
a single cyclic evolution, which means that the new phase is geometric because it depends on the winding 
number (if you go m times around a path, the phase is m times the phase for going once around the path). 

The main reason why people want to use geometrical phases for quantum computation is that frequently 
geometrical evolutions are easier to control and may also be more resistant to random noise coming from the 
environment. Since experimentally it is much easier to control the Hamiltonian than the actual state of the 
system, the adiabatic evolution is of importance. This means that the evolution of the system takes a long 
time compared to the characteristic dynamical time scale. This condition must be subtly changed in the 
framework of quantum open systems. The evolution must take a long time compared to the characteristic 
dynamical time scale but also it must be shorter than the decoherence time-scale introduced by the presence 
of the environment. In the model we have studied here, both conditions are fulfilled since decoherence is not 
effectively induced on the system for a realistic small environment. The results shown herein arc particularly 
relevant in the experimental realizations for the measurements of these phases. 
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